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Abstract 

The open quantum system formulation of the nuclear shell model, the so-called Gamow Shell 
Model (GSM), is a multi-configurational SM that employs a single-particle basis given by the 
Berggren ensemble consisting of Gamow states and the non-resonant continuum of scattering states. 
The GSM is of particular importance for weakly bound/unbound nuclear states where both many- 
body correlations and the coupling to decay channels are essential. In this context, we investigate 
the role of i=0 antibound (virtual) neutron single-particle states in the shell model description of 
loosely bound wave functions, such as the ground state wave function of a halo nucleus ^^Li. 

PACS numbers: 21.60.Cs,03.65.Nk,24.10.Cn,24.30.Gd 
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I. INTRODUCTION 



The theoretical description of strongly correlated open quantum systems (OQS), such as 
the weakly bound/unbound atomic nuclei or atomic clusters, requires the rigorous treatment 
of both the many-body correlations and the continuum of positive-energy scattering states 
and decay channels [30|. A major theoretical challenge is to consistently describe many- 
body states close to particle-emission thresholds where novel properties appear, such as, 
e.g., unusual radial features of halo states or threshold anomalies in the wave functions and 
associated observables. These features cannot be described in the closed quantum system 
(CQS) framework of a standard shell model (SM) which usually employs the single-particle 
(s.p.) basis of £^-functions of the harmonic oscillator. Representation of halo states in such 
a basis is not convenient. Moreover, the resonance and scattering states do not belong to 
the space of £^-functions. 

The solution of the configuration interaction problem in the presence of continuum states 
(the so-called Continuum Shell Model) has been advanced recently in the OQS formulation 
of the nuclear SM, such as the Shell Model Embedded in the Continuum 3j4| (real-energy 
continuum SM) and the Gamow Shell Model (GSM) S Q, 3 S H III, 12 , Q Q (SM 
in the complex fc-plane). A variant of the real-energy continuum SM, based on the Fesh- 
bach projection formalism and phenomenological continuum-coupling, has been proposed in 
Ref. Iia]. 



In the GSM. the multi-configurational SM is formulated|n the rigged Hilbert space Uflli: 
with a s.p. basis given by the Berggren ensemble [13, ll9[ consisting of Gamow (resonant 
or Siegert) states and the non-resonant continuum of scattering states. In general, resonant 
states correspond to the poles of the scattering matrix (S-matrix). These are the generalized 
eigenstates of the time-independent Schrodinger equation, which are regular at the origin 
and satisfy purely outgoing boundary conditions. The s.p. Berggren basis is generated by 
a finite-depth potential, and the many-body states are expanded in Slater determinants 
spanned by resonant and non-resonant s.p. basis states 0]. The configuration mixing 
induced by the GSM Hamiltonian assures a simultaneous treatment of continuum effects 
and inter-nucleon correlations. Consequently, the GSM, which is a natural generalization of 
the SM for OQS, gives a natural description of many-body loosely bound and resonant states. 
In general, among different poles of the one-body S-matrix, one takes into account bound 
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and decaying s.p. states to define the subset of resonant states in the Berggren ensemble. 
These poles, together with associated non-resonant continuum states, define the valence 
space. The actual selection of resonant states depends on the physical problem and on the 
convergence properties of resulting many-body GSM states. 

An important question concerning the GSM deals with the inclusion of antibound (vir- 
tual) states in the Berggren ensemble. Antibound states have real and negative energy 
eigenvalues that are located in the second Riemann sheet of the complex energy plane (the 



corresponding momentum lies on the negative imaginary axis) |20l . l2lL l22l |23[ . Contrary 
to bound states, the radial wave functions of virtual states increase exponentially at large 
distances. As often discussed in the literature, it is difficult to give a clear physical in- 
terpretation to virtual states. Strictly speaking, as the second energy sheet is considered 
unphysical and unaccessible through direct experiments, a virtual state is not a state but a 
feature of the system. If the virtual state has a sufficiently small energy, its presence has an 
appreciable influence on the behavior of the scattering cross section at low energies. Classic 
examples include the low-energy 15*0 nucleon-nucleon scattering characterized b y a large and 
negative scattering length , scattering of slow electrons on molecules 0, , and eep 



Coulomb system 



271 . Relatec 



to this is an increased localization of real-energy scattering 

states just above threshold 

Coming back to nuclear structure, it was argued that the neutron-unbound ^'^Li nucleus 
sustains a low-lying lsi/2 antibound state very close to the one-neutron (In) emission thresh- 
old j3] as a result of the inversion of Opi/2 and Is 1/2 shells j^ol. Although many theoretical 
calculations predict the Opi/2 — l'Si/2 shell inversion, the phenomenon still remains a mat- 



ter of debate 
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I . Experimenta. 
strength at the In-threshold in ^°Li 



ly, several 



33, 
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groups 
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have reported evidence of the £ = 



I ; however, no evidence of a weakly 
bound lsi/2 state has been found. Since, experimentally, the n+^Li has a large and negative 
scattering length, this may indicate the presence of the antibound lsi/2 state in ^°Li close 
to the In-threshold, though the presence of a low- lying Opi/2 state cannot be ruled out 13911 . 

nTi 

The GSM calculations that explicitly consider the antibound lsi/2 s.p. state I12| 
were performed recently for the g.s. of ^^Li. The authors argued that the presence of an 
antibound state was important for the formation of a neutron halo. They have also noted 
that the bound state wave function of ^^Li could be expanded in terms of the real-energy, non- 
resonant £=0 continuum, i.e., without explicit inclusion of the antibound state. Moreover, a 
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destructive interference between the lsi/2 antibound s.p. state and the associated complex- 
energy, non-resonant si/2 background was noticed. 

This work addresses the question of the £=0 virtual state for the description of a neutron 
halo by studying several complementary Berggren basis expansions and performing GSM 
calculations for ^^Li in different s.p. bases. In Sec. IHl the one-body Berggren complete- 
ness relations with and without the si/2 antibound state are compared in the i=0 channel. 
Section IIIll describes the g.s. of ^^Li in the schematic two-particle model using different 
s.p. Berggren ensembles, with and without the explicit inclusion of the antibound s.p. state 
lsi/2. The comparison of convergence properties of the corresponding GSM calculations 
with the number of Si/2 scattering states is made for the g.s. energy of ^^Li. A summary of 
results is given in Sec. IIVI 



II. COMPLETENESS RELATIONS INVOLVING ANTIBOUND si/2 STATES 



One-body Berggren completeness relations with bound and resonant states have already 
been extensively studied for neutrons 0, [l^ 



IS 



4G 







4l| and protons (7|,|42]. The standard 



Berggren completeness relation consists of a discrete sum over bound and resonant states, 
and an integral over non-resonant scattering states from the contour (see Fig. P): 



\un){un\+ / \u{k)) {u{k)\ dk 



(1) 



where a discrete sum runs over all bound (b) and resonant decaying (d) states lying above 
the complex contour . The continuous part takes into account the non-resonant scat- 
tering states lying on the contour. In the particular case of the £=0 neutron partial wave, 
there are no Si/2 resonances due to the absence of both Coulomb and centrifugal barriers. 
Consequently, a real-A; contour would have been sufficient to describe the si/2 neutron chan- 
nel. However, to investigate the convergence of imaginary part components of the expanded 
wave functions, we take the complex-fc contour close to the real-A; axis. 

In this work, we shall perform detailed studies of the Berggren expansion in a more 
general case when an £=0 antibound s.p. state is included in the basis. In this case, the 



Berggren completeness relation takes the form 111 l4ll|: 



\un){un\ + / \u{k)){u{k)\ dk 



(2) 



n£{a,b,d) 



Im[k] 




FIG. 1: Contours in the complex-/c plane used in the Berggren completeness relations for the Si/2 
(Sec. Oand Sec. and pi/2 (Sec. partial waves. The contour (OABC) is only used in 
the Si/2 channel; it allows wielding the antibound state lsi/2 (marked as 'a'). The contour 
(OA'B'C) is employed for Si/2 ^-^id •pxl2 channels and permits expansions of bound and resonant 
states only. 

where the sum also includes antibound (a) states lyin g ab ove the complex contour (see 



Fig. Q). It is to be noted that, as discussed in Ref. 4lJ, the contour is obtained by 
deforming continuously the contour placed in the fourth quadrant of the complex-fc 
plane so that it encompasses the antibound states of interest. In the unlikely situation that 
bound states of energies higher than antibound states (a) are present, they must be excluded 
from the sum 

Bound states can be expanded with either of the completeness relations (P) or Q , 
whereas antibound states can only be expressed through the more general completeness re- 
lation Indeed, only those poles of the 5'-matrix which are situated above the contour can 



be expanded using the Berggren basis 



4l| . The contour in the non-resonant continuum 



is usually discretized with a finite number of points in order to construct the Hamiltonian 

n 

matrix (see Ref. j6j for details). 
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TABLE I: Set of different WS potentials and lsi/2 states used in numerical tests of the Berggren 
completeness relation. The Berggren ensemble generated by a potential WS^''-' (second column) 
consists of the Osi/2 bound s.p. state, the contour in the non-resonant continuum (third column), 
and - possibly - the lsi/2 s.p. state (fourth column). 'No pole' denotes a situation where the 
virtual lsi/2 state is not included in the basis. In all the cases, the expansion has been carried out 
for the loosely bound lsi/2 s.p. state of ^^5*2. 



Case 


WS(o) 


Contour 


lsi/2 (WS(o)) 


(i) 


WSi 


Lt 


Antibound 


(ii) 


WSi 


Lt 


No pole 


(iii) 


WSs 


Lt 


Well bound 



In this study, the s.p. states entering the completeness relations and ((T)) have been 
generated using an auxiliary Woods-Saxon (WS) Hamiltonian: 

2m i + exp(^)' 

where m is the reduced mass of a neutron with respect to the ^Li core, (i=0.65fm, and 
i?o=2.7fm. The depth of the WS potential, Vo=50.5, 52.5, and 60.5 MeV, was adjusted 
to yield the lsi/2 eigenstate respectively antibound at —0.002955 MeV, loosely bound at 
—0.0329 MeV, and well bound at —1.0372 MeV. The corresponding WS potentials are 
denoted as WSi, WS2, and WS^ in the following. 

In order to test Berggren completeness relations (0) and (H)), we expand the lsi/2 eigen- 
state of a given WS potential in the basis generated by another WS potential (WS^^^^) of 
different depth: 

llsi/2) = ^c„|m„) + / c{k)\u{k)) dk (4) 

n J L+ 

where n is running either over bound and antibound \un) basis states in Q or over bound 
\un) basis states in (^, and L"*" is or L^, respectively. All the combinations of potentials 
studied are listed in Table HI The contours and used in this section are defined by 
vertices (all in fm"!): [O = (0.0, 0.0); A = (-0.01,0); 5 = (-0.01, -i0.02); C = (3.5,0.0)] 
and [O = (0.0, 0.0); A' = (0.1, -zO.Ol); B' = (1.5, 0.0); C = (3.5, 0.0)], respectively. 

To assess the quality of the Berggren expansion, we calculate the rms deviation from 
the exact lsi/2 halo wave function Uws2{^) WS2 obtained by a direct integration of the 



Schrodinger equation. The rms deviation is calculated separately for the real and imaginary 
parts, 



rms (Re [m]) 
rms(Im[M]) 



\ 



1 ^ 



i=l 



\ 



i=l 

for N=512 equidistant points on the real r-axis in the interval from r=0 to r=15 fm. In 
Eqs. dSI) and (jB)) Mws(o)(^) lsi/2 halo wave function of WS2 expanded in the basis 

WS(°) of Table E 

The calculated rms deviations (0) and (jH)) are shown in Fig. |21 One clearly sees that 
the Berggren basis containing an antibound state (case (i) in Table H} is less efficient in 
expanding a loosely bound lsi/2 state. The number of discretized scattering states in case 
(i) must be two-to-four times bigger than that in cases (ii) and (iii) in order to attain the 
same precision for the real part of the wave function. For the imaginary part, the difference 
is even more pronounced. 

Without an antibound state in the basis, 40 to 50 non-resonant scattering states are 
enough to obtain the precision of order 10~^ for the calculated si/2 wave function, whereas 
150 non-resonant scattering states are necessary to reach the same precision with this state 
included. In cases (ii) and (iii), one finds similar rms deviations, because the si/2 basis wave 
functions are in both cases either bound or close to the real (positive) A;-axis; hence their 
contributions add up constructively. On the contrary, in case (i), the virtual state with 
exponentially increasing wave function interferes destructively with non-resonant scattering 
states in order to produce the halo state. 

To illustrate this effect. Fig. El shows the distribution of squared amplitudes and c^{k) 
@. By construct,o„. the sum J] 4 + f c\m . always uormalized to one. As the 

contribution from the Osi/2 basis state is practically negligible, the set of reduces to the 
single real amplitude c\g^^^. In case (i), the non-resonant amplitudes c^{k) are much larger 
than their antibound counterpart c^^^^^, whereas in case (iii) the bound lsi/2 pole carries 
half of the halo wave function. Case (ii) lies between these two extremes, as the halo state 
is built entirely from the non-resonant continuum but in an essentially coherent way (cf. 
Refs. 




FIG. 2: Real and imaginary parts of the rms deviations ((3 for the Is 1/2 halo wave function 
expanded in different Berggren bases as a function of the number of ,31/2 scattering wave functions 
on L+. Cases (i)-(iii) of Tableware marked by dashed, solid, and dotted lines, respectively. 

The results displayed in Fig. |21 demonstrate that the inclusion of the antibound pole in the 
basis enormously enhances the role of the non-resonant continuum which has to efface the 
incorrect asymptotics of an antibound pole in order to create a bound state with decaying 
asymptotics. This behavior is opposite to what is found when including a narrow resonant 
state in the Berggren ensemble. Such a state always concentrates a fairly large part of the 
expanded wave function Q| . 
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FIG. 3: Real (circles) and imaginary (crosses) parts of the squared expansion amplitudes and 
c^{k) @ corresponding to the Berggren bases (i-iii) of Table in the complex A;-plane for the 
lsi/2 halo wave function. The amplitude Cqsi/2 practically negligible. In cases (i) and (iii), the 
amplitude cf^^^^ (magnified by a factor of 10) is represented by an arrow placed arbitrarily at A; = 
0. 

III. CONTRIBUTION OF ANTIBOUND STATES TO THE TWO-BODY HALO 
STATES: EXAMPLE OF ^^Li 

In order to assess the influence of antibound states on many-body halo states, we shall 



111 ll^ of two valence particles, coupled to J'^=0+, moving 
outside the inert core. In particular, we shall view the g.s. of ^^Li in terms of two neutrons 



employ a schematic model 



in s and p angular momentum states coupled to the ^Li core. Our aim is not to exactly 
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reproduce the structure of ^^Li but to learn how antibound s.p. states enter a two-body halo 
wave function. 

The nuclear Hamiltonian is given by a WS potential representing a ^Li core, to which a 

n 

Surface Gaussian Interaction (SGI) [7] is added, modeling the residual interaction between 
the two valence nucleons: 

^2 



H 



i 

■ exp 



£-v„./(n)-4v.„(5-.-:)^ 



<i/(r.) 



dr 



i<j 



■5(|rl| + |r;-|-2i?sGi). 



(7) 



Here, we employ the set of WS parameters WSi of Sec. HH and ^^0=21.915 MeV is the 
strength of the spin-orbit potential. For these parameters, the lsi/2 state is antibound and 
the Opi/o state is a resonance with the energy £'=0.24 MeV and the width r=118 keV, in 

n n 

fair agreement with experimental data |29|, |39| . The range of the residual SGI interaction is 
/i=l fm, i?sGi=4.4 fm, and Vsgi=— 1196 MeV fm^. These parameters have been adjusted 
to reproduce the binding energy of ^^Li with respect to ^Li (Eb =-0.295 MeV), and a nearly 
50% weight of the component, as suggested by the data 29|. 

The valence space consists of si/2, P3/2, and pi/2 neutron wave functions, from which the 
Osi/2 and 0^3/2 bound states have been removed as they are part of the ^Li core. The P3/2 
s.p. space contains the 0^3/2 bound state pole and the real-energy ^3/2 contour with k^ax='^ 
fm^^. The pi/2 space consists of the Opi/2 resonant pole and the pi/2 contour of an type 
(see Fig. d) defined by the vertices (all in fm"!): [O = (0.0, 0.0); A' = (0.1, -^0.1); = 
(1.0, 0.0); C = (2.0, 0.0)]. The Si/2 space contains either a scattering component (L^) or the 
lsi/2 antibound pole to which the si/2 contour of an L+-type is added. The contours 
and in the Si/2 channel are defined by the vertices (all in fm^^): [O = (0.0, 0.0); ^4 = 
(-0.01, 0);B = (-0.01, -i0.04); C = (2.0, 0.0)] and [O = (0.0, 0.0); A' = (0.1, -iO.Ol); B' = 
(1.0, 0.0); C = (2.0,0.0)], respectively. The P3/2 and pi/2 contours are discretized with 30 
and 32 points, respectively. Each point represents one shell in GSM calculations. For this 
level of contour discretization and the momentum cut-off, the theoretical error on energies 
is about 1 keV for the real part and 0.01 keV for the imaginary part. 

The convergence of the ^^Li g.s. energy as a function of the number of Si/2 non-resonant 
scattering shells is shown in Table |nj The number of Si/2 shells on each segment of the 
contours and is chosen so as to minimize a spurious width of the g.s. It is seen that 
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TABLE II: Energy E (in MeV) and width T (in keV) of the ^^Li g.s. function of the number 
Ng^^^ of the non-resonant scattering shells in the discretized Sx/2 continuum. The values of E/T[L'^] 
(second and fourth columns) are calculated with the Sx/2 space consisting of the isi/2 antibound 
state and the associated si/2 contour. The values of E/T[L'^] (third and fifth columns) are 
obtained with the Si/2 space consisting of the Si/2 contour only. See Fig. ^ for the definition 
of different contours. 





E[L+] 








10 


-0.314 


-0.291 


65.274 


-3.644 


20 


-0.292 


-0.295 


2.307 


0.025 


30 


-0.294 


-0.295 


0.876 


-0.003 


40 


-0.294 


-0.295 


-0.425 


-0.007 


50 


-0.295 


-0.295 


0.075 


-0.009 


60 


-0.295 


-0.295 


-0.005 


-0.009 



the energy converges much faster when the antibound state is not present in the s.p. basis. 
In this case, the calculated g.s. energy of ^^Li attains a precision of 0.1 keV with 20 shells 
included, whereas as many as 50 shells are necessary to obtain the same level of precision if 
the antibound pole is present. 

It is instructive to inspect the g.s. wave function of ^^Li expressed in different Berggren 
bases. Table lllll compares the GSM results obtained in the Berggren basis including the 
lsi/2 antibound pole and Si/2 non-resonant scattering states from the contour (case (i) 
of Table lU with those obtained using the contour for the Si/2 part (case (ii)). (The 
P3/2 and pi/2 spaces are defined as above.) To determine precisely the valence neutron 
configurations, 60 shells were taken along each of the scattering contours: Si/2, Pi/2, and 
^3/2- The expansion amplitudes of the p components are identical in both cases and tote 
up to about 49%. As expected, important cancellations appear in case (i). Here, about 
50% of the sum of squared amplitudes for configurations with two neutrons in Si/2 non- 
resonant scattering states is cancelled by contributions from configurations involving one 
neutron in a scattering state and another one in the lsi/2 antibound state. The antibound- 
antibound configuration (lsi/2)^ contributes only ~10%. In case (ii), however, all scattering 
two-neutron configurations add up coherently. 
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TABLE III: Real and imaginary parts of squared amplitudes of the GSM configurations involving 
neutrons in si/2 orbits in the g.s. of ^^Li. The GSM calculations were performed using two Berggren 
bases: (a) including the lsi/2 antibound pole state and the non-resonant scattering states Si/2 along 

(second and fourth columns), and (b) including the non-resonant scattering states S1/2 along 

(third and fifth columns). For the definition of contours, see Fig. ^ 



Configuration 


Re[c2]{L+} 


Re[c2]{L+} 


Im[c2]{L+} 


Im[c2]{L+} 


(1^1/2)' 


0.0990 




-9.6033-10-6 




(lSl/2 S1/2) 


-0.5887 




2.3369-10-5 




(^1/2)' 


1.0034 


0.5137 


-8.0720-10-6 


-1.4650-10^5 



IV. CONCLUSIONS 

The GSM has proven to be a powerful approach for the microscopic description of loosely 
bound and resonant states. The fact that the underlying Berggren basis directly incorpo- 
rates one-particle continuum and the proper treatment of many-body correlations through 
configuration mixing makes it a perfect tool for a theoretical description of loosely bound 
many-body states such as nuclear halos. 

In this study, we investigated the importance of including the virtual i=0 state in the 
s.p. basis for a modeling of one-body and two-body neutron halos. This question has 
experimental relevance: the data seem to suggest that the presence of the si/2 antibound 
state in ^°Li is correlated with the appearance of the two-neutron halo in ^^Li. 

Our calculation suggests that neither conceptual nor numerical gain is achieved if the 
antibound states are directly considered in the Berggren ensemble. In both one- and two- 
particle cases, it is always more advantageous to use a s.p. basis which contains non-resonant 
scattering states and, possibly, a bound S1/2 pole. As the halo wave function has a decaying 
character at large distances due to the exponentially increasing asymptotics of the virtual 
state, its inclusion in the basis always induces strong negative interference with the scattering 
states. Therefore, adding antibound states to the basis is not beneficial, as more discretized 
scattering si/2 states are necessary to reach a required precision without providing any new 
information about the many-body wave function. 
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The use of antibound states in a Berggren basis can only be justified if the many-body 
state has a virtual character. This is certainly not the case for bound states, including 
halos. For instance, the g.s. halo in ^^Li is well described by using the Berggren basis solely 
involving the S1/2 non-resonant scattering continuum and choosing the integration contour 
close to the real-A; axis. While the presence of an antibound state does require an increased 
density of discretized states around k=0, this can be handled efficiently by employing the 
density matrix renormalization group 

better numerical precision by suppressing the cancellations and leaving all the physical 
properties unchanged. Hence, the use of antibound states in GSM becomes questionable, as 
it worsens spurious effects due to continuum discretization. 
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